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Abstract: The self-couplings of the electroweak gauge bosons are probed at hadron col- 
liders through the production of a massive gauge boson and a photon. To extend the 
theoretical description of this type of final states towards next-to-next-to-leading order 
(NNLO) in QCD, we derive the two-loop QCD corrections to the helicity amplitudes de- 
scribing the production of a massive gauge boson in association with a real photon. Our 
results are obtained by applying projectors to the general parton-level tensor structure. 
The leptonic decay of the vector boson is included, thus allowing for a fully exclusive de- 
scription of the final state. The infrared poles of the amplitudes are described by an infrared 
factorization formula. We provide an analytic expression for the finite remainder of the 
amplitude in terms of one- and two-dimensional harmonic polylogarithms. The amplitudes 
are expressed in the physical kinematics relevant to gauge-boson-plus-photon production 
at hadron colliders. As a by-product, we also derive the two-loop QCD amplitudes for 
vector-boson-plus-jet production at hadron colliders. 
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1. Introduction 

Pair production of electroweak gauge bosons (7, W ± , Z°) offers a wide spectrum of observ- 
ables, which allow to test the theory of the electroweak interaction, to probe the Higgs 
mechanism of electroweak symmetry breaking and to search for physics beyond the stan- 
dard model. While photons are directly observed in the detector, the massive W and Z 
bosons are identified from their leptonic decay modes. 

The standard model predicts specific values and structures for the couplings among the 
electroweak gauge bosons: W ± , Z° and 7. Physics effects beyond the standard model could 
modify these gauge boson self-couplings [1,2]. Observations of such anomalous couplings 
may help to constrain new theory models and could provide indirect evidence for new 
physics effects at energy scales above the nominal collision energy. The couplings of the 
massive , Z® bosons to the photon are determined at hadron colliders by measuring 
W ±n f and Z°7 production cross sections and comparing them to theoretical predictions. 
Measurements have been carried out at Tevatron [3], and first results from the LHC are 
already becoming available [4]. The accuracy of the coupling determination is potentially 
limited by both the experimental accuracy and by uncertainties inherent to the theoretical 
prediction. 

At present, W^j and Z°~/ production at hadron colliders is described theoretically to 
next-to-leading order (NLO) in QCD [2,5,6] and to NLO in the electroweak theory [7]. With 
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increasing order in perturbative QCD, new production channels (with new combinations 
of parton distributions) for vector boson pairs start contributing; the complete spectrum 
of partonic channels is only present from next-to-next-to-leading order (NNLO) onwards. 
Moreover, the inclusion of NNLO corrections to gauge boson pair production will lower 
the inherent theoretical uncertainty of the prediction (usually quantified by variation of 
the renormalisation and factorisation scales) and allow for a fully consistent inclusion of 
NNLO parton distribution functions. 

The calculation of gauge boson pair production at NNLO requires three types of ingre- 
dients: the two-loop partonic 2 — > 2 matrix elements for the production of the gauge boson 
pair under consideration, the one-loop partonic 2 — > 3 matrix elements for the production 
of the gauge boson pair in association with an extra parton and the tree-level 2—7-4 matrix 
elements involving two extra partons. The latter two contributions are equally contributing 
to the NLO corrections for the production of a vector boson pair with an extra jet, which 
have been computed for 77J [8], Vjj [9] and VVj [10] already some time ago. At NNLO, 
the contributions from both these channels will contain infrared singularities from one or 
two final state partons becoming soft or collinear. These singularities cancel only when 
combined with the infrared-singular two-loop contributions, such that a method is needed 
for their extraction from the real radiation processes. Several methods have been applied 
successfully in NNLO calculations of exclusive observables in the recent past: sector de- 
composition [11], (/^-subtraction [12] and antenna subtraction [13-15]. It should be noted 
that the (/T-subtraction method is restricted to observables that are described by non-QCD 
processes at leading order, which is the case for vector boson pair production. The first 
calculation of NNLO corrections to a vector boson pair production process (pp —> 77) has 
been performed most recently [16] using this method. 

The two-loop parton-level matrix elements for vector boson pair production are at 
present known only for qq — > 77 [17] and gg — > 77 [18] (where the latter two-loop amplitude 
formally contributes only beyond NNLO). The high energy approximation for qq — > W + W~ 
and qq — > Z°Z° has also been derived [19]. It is the purpose of this paper to compute the 
two-loop corrections to the matrix elements for the production of a massive vector boson 
and a photon: qq — > W^j and qq — > Z°7. The calculation follows closely the techniques 
that were employed in the calculation of two-loop corrections to the 7* — > qqg matrix 
elements [20,21], which were a crucial ingredient to the NNLO corrections [22,23] to three- 
jet production and related event shapes at e + e~ colliders. 

This paper is structured as follows: in Section 2, we fix the notation and discuss the 
basic helicity structure of the process under consideration. The calculation of the two-loop 
amplitudes is described in Section 3, and the results are discussed in Section 4. The two- 
loop helicity amplitudes are obtained in a closed analytic form, which is however too large 
to be quoted in the paper, and we enclose computer algebra files containing the results 
with the submission. We performed several non-trivial checks on the results, which are 
described in Section 5. We conclude with an outlook in Section 6. We enclose appendices 
with the one- loop helicity amplitudes, examples of selected colour coefficients in the two- 
loop amplitudes and with a discussion on kinematical crossings, including the two-loop 
matrix elements relevant to vector-boson-plus-jet production: qq — > Vg and qg — > Vq. 
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2. Kinematics and basic helicity structure 



The production of a massive vector boson and a photon in quark-antiquark annihilation 
is related by crossing to the decay of the vector boson into a quark-antiquark-photon final 
state, which has the same kinematics as three-jet-production (3j) in e + e _ annihilation. 
Technically, the calculation of QCD corrections to the qq — > V7 amplitudes is thus similar 
to previous calculations for the helicity amplitudes for 3j-production, which have been 
derived to two-loop accuracy in QCD [21]. 

Including the leptonic decay of the vector boson, the partonic subprocesses yielding 
Vj final states are: 

q(P2) + q(pi) -> 7(-P3) + Z°(q) -> 7 (-p 3 ) + Z+(p 5 ) + r (pe) , 
?(P2) + q'ipi) -»• 7(-Ps) + ->• 7(-Pa) + "(ps) + i~(pe) , 

?(P2) + ^(Pi) -> 7(-Ps) + -»• 7(-p 3 ) + l + (P5) + v(pe) • 

The Z°-boson process implicitly includes also a contribution from an off-shell photon 7*. 
In the most general case where two quarks of two different flavours appear, to fix the 
conventions, we will refer from now on to p± as the momentum of the anti-quark q' and as 
P2 to the momentum of the quark q. The momentum of the vector boson is given by 

9"=l£+pff. (2-1) 
It is convenient to define the invariants 

«12 = (Pi +P2f , S13 = (Pi +P3) 2 , S 2 3 = (P2 +P3) 2 , (2-2) 

which fulfil 

Q 2 = (Pi + P2 + P3) 2 = S12 + S13 + S23 = S123 , (2.3) 
as well as the dimensionless invariants 

^ = S12/Sl23, y = Sl3/si23, ^ = S23/si23, (2-4) 

which satisfy x + y + z = 1. 

In 3j-production, is time-like (hence positive) and all the Sij are also positive, which 
implies that x, y, z all lie in the interval [0; 1], with the above constraint x + y + z = 1. For 
the V7 production, q 2 remains time-like, but only s±2 is positive: 

q 2 > , S12 > , S13 < 0, S23 < , (2.5) 

or, equivalently, 

x>0, y<0, z<0. (2.6) 

It was shown in [24] that the kinematical situation of this configuration can be expressed 
by introducing new dimensionless variables 

S13 y q 2 1 

u = = , v = — = -, 2.7 

S12 X S12 x 
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which fulfil 



< u < v. 



0<v<l 



The helicity amplitudes for qq — > can be expressed as a product of a partonic 
current and a leptonic current L M : 



A(P5,Pe;Pi,P3,P2) = L^(p 5 ;p 6 )Sf J ,{p 1 ;p3;p 2 ) 



(2.8) 



Only the partonic current receives contributions from QCD radiative corrections, and it 
can be perturbatively decomposed as: 



Sa(pi;P3-,P2) = V^TTO 



Pi;P3-,P2) + (7^) Sjp(p 1 ;p 3 ;p 2 ) 
+ (^ys^( Pl ;p 3 ;P2) + 0(a 3 s )). 



(2.9) 



It is a colour-singlet. The vector boson decay to a lepton-antilepton pair is described by 
a leptonic current. To be as general as possible, we consider only the basic amplitude 
structure in the partonic and leptonic current, and include charges and coupling factors 
related to the massive vector boson only when assembling the final results. We have 
extracted a factor e = \/Aira for the photon coupling in the partonic current, such that all 
the quark charges will be expressed in units of e. 

The most general structure of the partonic current can be derived from symmetry 
considerations [21]: 



S»(pi;p3;p2) 



An Tn^ + A12 Ti2 M 



+ ^13 ?W 



+ A2I T2ln + A22 T22/M + A23 T23 /1 



(2.10) 



where and T M are the following tensor structures: 

Tijv, = v(pi) fa e 3 • pi p iM - p jfl + Sj ^hhl^ u(p 2 ), 



v(pi) 
T„ =v{pi) 



h £3 • P2 Pjfi 



2 

S23 



2 r 3Pj l x+ Sj £lM3 
1 



u(p 2 ), 



S23 1^3 ■ Pi + T^hhl^ - s 13 1^3 ■ P2 + ^l^hh 



u(p 2 ), (2.11) 



where we defined: 



SU = Si2 + Si 3 , S 2 4 = Si2 + S23, «34 = «13 + «23- 

The tensor coefficients Aij and B can be determined by appropriate projectors, applied 
to the Feynman-diagrammatic expression of the amplitude. Projections on the diagrams 
are performed in dimensional regularisation in d = 4 — 2e dimensions. The projectors can 
be found in [21]. 

Each of the unrenormalised coefficients and B has a perturbative expansion of the form 



4 un 

A i3 



\/4~ 
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ira 



B (0),un + f ««\ B (l),un + ( ^±\ 2 B (2),nn + ^3) 

\2ttJ \2tt J 



(2.12) 



where the dependence on (513,523,5123) is implicit. 

By fixing the helicities of the partons, the partonic current can be cast in the usual 
spinor helicity notation [25]. All helicity configurations can be obtained from the amplitude 

S^PtPt) = ^ < 12 >[ 13 ] 2 (Pi^n +P2,A l2 + p 3 ,A 13 ) - ITmI 2 ) S ^B 



+ ^=[13] [3 \^\ 2) 



S2$B + - {{A u + A 12 )s 12 + {A u + i4i 3 )si 3 + (A i2 + Ai 3 )s 2 3) 



(2.13) 



by charge and parity conversion. For q'(pi), q{p 2 ) incoming, the above amplitude corre- 
sponds to a right-handed current. Notice that (2.13) has been obtained assuming that the 
momentum of the photon is — p 3 . We have: 



Sl(pi;p 3 ;p 2 ) = [Sr{pi;p 3 ;p 2 )]* > 
sKpt\pt ;p 2 ) = -SR{p 2 ;pt;pt), 
Sr{pi\p 3 ;p 2 )= [Sr(p 2 ;p 3 ;Pi)Y 



(2.14) 



It is also straightforward to include the spin-correlations with the leptonic decay products 
by contracting the partonic current with the leptonic current for fixed helicities of the 
final state leptons. Consider the decay of the vector boson V into two leptons: 

V(q)—H + (p 5 ) + r(p6). 

The purely vectorial tree- level leptonic current reads: 

L»(p 5 ,IHi) = u( P e)-fv( P5 ), (2.15) 

where in the case of an outgoing lepton-antilepton pair L fl {p^ ,p§) corresponds to a right 
handed current, and L^p^p^) to a left-handed current. We find straightforwardly: 



L>* R {Ps,pt) = [6 |7l 5), Ll{pi,p 6 ) = [5 | 7 1 6) = [L^, %+)]*. 



(2.16) 



In order to write down the lepton-parton contraction it is convenient to introduce the 
set of helicity coefficients defined in [21]: 



a{u, v) = 
(3{u,v) = 
j(u,v) = 



si 3 s 23 



£13 
4 

S13S23 



2S + Ai2-An), 
2S23-B + 2{s 12 + S13M11 + s 23 {A 12 + A 13 fj . 
A u - A13), 



d{u,v) = — An, 



(2.17) 
(2.18) 
(2.19) 
(2.20) 
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which, from their definition in terms of the coefficients A^ and B, respect the relation 

2S123 



a(u, v) — f3(u, v) — j(u, v) 



Sl2 



-S(u, v) 



0. 



(2.21) 



The relations above can be inverted for An, A\ 2 + 2B and ^13, and in these variables the 
contracted amplitude assumes a particularly simple form. We take the contraction of the 
right-handed quark current with positive photon helicity, and the right-handed leptonic 
current, as basic object from which all other helicity configurations are obtained: 



a rr(P5,Pg;Pi,P3,P2) = L R {p 5 ;p£ )S R:IJi (p 1 ;pf;P2) 



= -2V2 



<25)<12)[16] 
(13) (23) 



a(u, v) 



(25) [36] (15) [13] [36] 



(13) 



(13) [23] 



(2.22) 



The unrenormalised helicity amplitude coefficients a, (3 and 7 are vectors in colour 
space and have perturbative expansions: 



\2ir/ \2ir J 



(2.23) 



for ft = a, /3, 7. The dependence on (u, v) is again implicit. 

From A~ft R (p5,pQ;pi,p3,p2), all other helicity amplitudes can be obtained by parity 
and charge conjugation. Axial contributions from the weak gauge boson couplings can be 
accounted for in a straightforward manner, by simply reweighting the different right-handed 
and left-handed helicity amplitudes with appropriate weights. 

The eight possible helicity configurations are obtained from A^ R as follows: 



l r (ps ;pt)SR»(pi;pt;P2) 

L R (P5;pt)S R ^;p 3 ;p^) 
L R (P5;pt)S Lll (pi;pp,p 2 ) 
L R (P5;pt)S Lf i(pi;p 3 ;P2) 



A RR (p5,Pa;Pi,P3,P2), 
A RR (p5,p 6 ;pi,P3,P2) -- 

A RL (p 5 ,P6-,Pl,P3,P2) = 
A R l(P5iP6'>P1iP3iP2) = 



[- A RR(Pe,P5-,P2,P3,Pl)]*, 
- A Rr(.P5,P6-,P2,P3,Pi), 

[ A rr(P6,P5;Pi,P3,P2)]*, (2.24) 



A t R (P5,P6',PUP3,P2) = A RR (P6,P5-,Pl,P3,P2), 
A l R (P5,P6-,Pl,P3,P2) = [- A RR (P5,P6-,P2,P3,Pl)]*, 
A l L (P5,P6;Pl,P3,P2) = - A RR {P6,P5-,P2,P3,Pl), 
A I L (P5iP6;PliP3,P2) = [ A RR (P5,P6;Pl,P3,P2)]* ■ (2-25) 

The general form of the gauge boson coupling to fermions is: 

v v,/i/ 2 = _ ieV VJih with e = ^/4^ ; (2.26) 

whose explicit form depends on the gauge boson, on the type of fermions, and on their 
helicities: 

K' h " = > (i^) + Rj, h t,. (±±*) • (2.27) 



L Z(pf-,P6 )SR»(Pi;pt;pt) = 

L l l(pi;PQ)S R ^;p 3 ;p^) = 
L /1 L (pt;Ps)S L ^;pp 1 p 2 ) = 

Ll(pt;P6)S L ^{pi;P3;P2) = 
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The left- and right-handed couplings are identical for a pure vector interaction, and are 
in general different if vector and axial-vector interactions contribute. Their values for a 
photon are 

L lfr = R lh = - e h S hh' ( 2 - 28 ) 

while for a Z boson 

z _ lj l - sin 2 6 w e h z smO w e fl 

L h h ~ „:„/. — d hh > R h h - — iits — °hh ' I 2 - 29 ) 



and finally for a W ± 

1 



L ^ = VT^F' fl "= a (2 ' 30) 



The charges ej are measured in units of the fundamental electric charge e > 0. 
The vector boson propagator can be written as: 



p ^ ) = ! i||f' (2 - 31) 



->v i„ c \ — > xv 



where Ay„(g, £) and Dy(q) are, respectively, the numerator and the denominator in the 



R$ gauge: 

= (-9»„ + (1 - O^^t) , (2-32) 

D Z ,W± (<?) = {q 2 -M 2 + iT v M v ) , (2.33) 
D^q) = q 2 . (2.34) 

In the narrow-width approximation we can simplify expression (2.33) to 

D z ,w±(<l) ~^VMy and q 2 = My, (2.35) 

where My is the mass of the vector boson, while Ty is its decay width. 

Since we do not consider any electroweak corrections, the vector boson V is always 
coupled to a fermion line and this allows us to neglect the dependence (or equivalently 
to put £ = 1). With this notation we obtain for the different choices of V = (j*,Z, H /± ), 
and helicity combinations (with the obvious notation pij = Pi+Pj): 

R v R v 

Mv(Ps,Pt->Pl>Pt>Pt) = -*(47ra) q ' q \ f f A+ R {p 5 ,p 6 ;p 1 ,p 3 ,p 2 ), (2.36) 

jV R V 

My(ps,pt;pt,P3,P2) = -^(4vra) q ' q \ f f [A+ R (p 6 ,p 5 ; Pl ,p 3 ,p 2 )}* , (2.37) 

jV ftV 

M v (p5,pi;pi,p+,p 2 ) = -i(4vra) q ' q ' 2 f f [-A+ R (p 5 ,p 6 ;p 2 ,p 3 ,p 1 )}, (2.38) 

R v R v 

Mv(j>5,pt;pT,I%,Pt) = ~H^a) q ' q2 , f f [-A+ R (p 6 ,p 5 ;p 2 ,p 3 , Pl )]*. (2.39) 

The corresponding amplitudes for left-handed leptonic current can be obtained simply 
interchanging ^5 o p§ and -Rj r j 6 — > -^j%/ 6 - 
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^WW v 




(2) 



V 




(3) 



Figure 1: Abelian and non-abelian tree- level contribution. 



3. Outline of the calculation 

The two-loop corrections to the coefficients Q for W^j and Z°j production can be eval- 
uated through a calculation of the relevant Feynman diagrams. The diagrams which con- 
tribute at the two-loop level can be organised in different classes, some of which are present 
in both V = Z°, cases, while others contribute only to one of the two processes. 

To identify the different classes of diagrams, it is useful to start with the tree-level. 
We can define three classes of processes, each represented by a single diagram: 

1. We call the contribution from diagram (1) in Figure 1, where the photon is 
attached on the quark q'. The charge factor of this diagram is e q i. 

2. We refer as to the contribution from diagram (2) in Figure 1, where the photon 
is attached on the quark q. The charge factor of this diagram is e q . 

3. ^3 is finally the contribution from diagram (3) in Figure 1, where an off-shell W* 
radiates the final state. The charge factor of this diagram is unity. 

The &i are measured in units of e, so that e q > — e q = 1. 

One can compute the through the use of the projectors defined above, and once 
the three different contributions are known one can reconstruct the correct values for the 
helicity coefficients as: 



n 



(o) 



U, 



qq' 



u, 



qq' 



n\ 0) + or ) + 



(0) 



>(0) 



n 



(0) 



(3.1) 







(o) 



e a ( n\ 0) + n 



,(0) 



(3.2) 
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where Uij are the CKM matrix elements. To simplify the expression above, we made use 
of the fact that the Z° boson does not couple to the photon, and that in this case there is 
no flavour change, i.e. e q i = e q . 

At leading order, the d-dimensional expressions for the coefficients read: 



4(d - 3) V u uj ' 2 4(d - 3) V u u 

col d — 4 / v 1 \ it , 
- '2 + • (3-3) 



4(d-3) V « «/ 1-w 



73 



(o) _ d - 4 



(I + 3 . (3,, 



4(d - 3) V u 
Using (3.1) and (3.2) we find 

= U qq < (e q , - ^) , C = U qq , [e q , - ^) , Z ] = 0, (3-6) 

4 0) r = e «' 4%=^, 7^=0. (3.7) 

At one-loop, the same classification of contributions applies: , A further 

type of diagrams, with both the photon and the gauge boson which couple to a closed 

(2) (2) (2) 

quark loop, is zero due to colour conservation. At two loops, besides I\ ,1^ ,1% , three 
further classes of diagrams appear: 

(2) 

1. I\ are the diagrams where both 7 and V couple to the same fermion loop, as depicted 
for example in Figure 2, diagram (4). 

This contribution is denoted by Npy and is proportional to the charge weighted sum 
of the quark flavours. In case of a 7* exchange we find 

V p 2 

N F „ = . (3.8) 

Considering Z- interactions, the same class of diagrams yields not only a contribution 
from the vector component of the Z, which for the right-handed quark amplitude is 
given by 

Y.q> ( L q >gi + R q>q>) e q> 
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TOOT 

(4) 



(VW 





(6) 




Figure 2: Examples of two- loop diagrams in the classes I\ , Jg and 1$ 



but also a contribution involving the axial couplings of the Z. This contribution van- 
ishes identically for Z°7 production, already before summing over the quark flavours 
inside the loop. In the case of exchange charge conservation ensures that 

N F>W ± = 0. (3.10) 

(2) 

2. ig are the diagrams where the photon couples alone to a fermion loop, while V 
couples to the fermion line, as depicted in Figure 2, diagram (5). This class of 
diagrams has to sum to zero due to Furry's theorem. 

(2) 

3. Iq are finally the diagrams where V couples alone to a fermion loop, while 7 couples 
to the fermion line, as depicted in Figure 2, diagram (6). These diagrams give both 
a vector and an axial contribution, where the vector contribution is again zero due 
to Furry's theorem, while the axial contribution is zero in the framework of massless 
QCD [21]. 

(2) (2) 

We explicitly evaluated the contributions from classes Ik , Jg , and their vanishing 

provides a check on our calculation. 

(2) 

The classes ij , j = 1,6 exhaust all the possible two- loop QCD diagrams which can 
contribute to the production of a pair V7, whatever is the identity of the vector boson V. 
When computing the helicity amplitudes, one can evaluate the contributions from these 
six classes of diagrams independently, without keeping track of the axial current contribu- 
tions, thus considering the vector boson V as an off-shell purely vector particle. Once the 
(2) 

I- are known, as for the tree-level case, one can reconstruct the proper amplitudes for 
V = Z°, summing these six contributions up, multiplied by appropriate weights. 
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The calculation proceeds as follows. The 143 two-loop diagrams belonging to the 
classes j = (1,2,4,5,6) are produced using QGRAF [26] while, on the other hand, we 

(2) 

did not need to evaluate explicitly the diagrams in class 1% , since they only account for 
the QCD corrections of the quark form factor, which is known up to three loops in the 
literature [27,28]. 

The tensor coefficients are then evaluated analytically diagram by diagram applying 
the projectors defined in [21]. As a result, one obtains the tensor coefficients in terms of 
thousands of planar and non-planar two- loop scalar integrals, which can be easily classified 
in two auxiliary topologies, one planar and the other non-planar [20]. Through the usual 
IBP identities [29] one can reduce independently all the integrals belonging to these two 
auxiliary topologies to a small set of master integrals. This reduction is performed using 
the Laporta algorithm [30], implemented in the Reduze code [31]. All the masters for such 
topologies are known [32] as series in the parameter e = (4 — d)/2, through a systematic 
approach based on the differential equation method [32,33]. The masters are expressed as 
Laurent expansion in e, with coefficients containing harmonic polylogarithms (HPLs, [34]) 
and two-dimensional harmonic polylogarithms (2dHPLs, [32]). Numerical implementations 
of these functions are available [35]. For all the intermediate algebraic manipulations we 
have made extensive use of FORM [36] . 

The two-loop unrenormalised helicity coefficients Q^' nn can then be evaluated as linear 
combination of the tensor coefficients, and in particular they can be evaluated separately 
for every class of diagrams: 



As for the tree level case, it is trivial to reconstruct the amplitudes for the processes 
considered as linear combinations of these six amplitudes. 

We start considering the case where V = W^ 1 . The W boson couples only to left- 
handed fermions and charge conservation implies that Npw = 0- The amplitudes for 
VF ± 7 production at two loops thus receive contributions only from three of the six classes 
of diagrams above, i.e.: 



As for the tree level, Uij are the CKM matrix elements. 

For the V = Z°,j* case, = at all orders, since the Z° and 7* are electrically 

(2) 

neutral, while I\ is non- vanishing since charge conservation does not forbid such diagrams 
anymore. 




(2),un 



with j = 1, 6 . 



(3.11) 




(3.12) 




(3.13) 



for V = (Z°, 7 *). 



4. Two-loop helicity amplitudes 



Renormalisation of ultraviolet divergences is performed in the MS scheme by replacing the 
bare coupling ao with the renormalised coupling a s = a s (/U 2 ), evaluated at the renormali- 
sation scale /U 2 . Since the tree amplitudes are of 0(oP s ), we only need the one loop relation 
between the bare and renormalised couplings: 



Q /U b e = a s /i 



'-T(£)+°e» 



(4.1) 



where 



S e = (4vr)V 



■«7 



with Euler constant 7 = 0.5772 . . . 



and /Xq is the mass parameter introduced in dimensional regularisation to maintain a di- 
mensionless coupling in the bare QCD Lagrangian density. /3q is the first coefficient of the 
QCD /3-function: 



6 



with the QCD colour factors 

C A = N, 



c F 



N 2 



2N 



T R = -. 



(4.3) 



The renormalisation is performed at fixed scale /U 2 = q 2 . The renormalised helicity coeffi- 
cients read: 



n (0) = n (0),un ) 

n {1) = sr 1 ^'™, 



(4.4) 



The full scale dependence of the coefficients can be recovered from the renormalisation 
group. It reads: 



Q.(n 2 ,a s (n 2 )) = V&raSij 
'a s (n 2 



a s (p 2 ) 
2ir 



+ 



2tt 



n {2) + (3 n {1) in 



/i- 



+ 0{a\) 



(4.5) 



After performing ultraviolet renormalisation, the amplitudes still contain singularities, 
which are of infrared origin and will be analytically cancelled by those occurring in radiative 
processes of the same order. Catani [37] has shown how to organise the infrared pole 
structure of the one- and two-loop contributions renormalised in the MS-scheme in terms 
of the tree and renormalised one-loop amplitudes. The same procedure applies to the 
tensor coefficients. Their pole structure can be separated off as follows: 



ft« = /( 1 )(e)0(°)+0( 1 )- finite , 
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= ^-Ij(D( e )/W( e ) _ ^jw (e ) + e -7£0__^) ^ + ^ jd) (2e) + H (2 )(e) j n( o) 

+I«(e)^ 1 )+ft( 2 )' finite , 



where the constant i^T is 



K 



67 vr 2 \ 10 

is -y c *-y TrNf - 



(4.6) 



(4.7) 



In our case, there is only a quark-antiquark pair present in the initial state, so that 
is given by, 



/ (1) ( £ ) = 



=^7 



2r(l - e) 



N 2 — I ( 2 3 



(4.8) 



where, since we have set ji 2 = s 2 23 : 



'12 



2 \ e 



Sl2 y 



(4.9) 



Note that on expanding S12, imaginary parts are generated, the sign of which is fixed by 
the small imaginary part +i0 of s\2- The hard radiation constant is a scalar in colour 
space: 



4eT(l - e)' 



with 

where in the MS scheme 

Hp = (N 2 - 1) 



7 409 llvr 2 \ 2 /3 3 ir 2 
4 C3 + 864 ~~96 _ J + \ 2 32 ~ ~8 J N 2 



(4.10) 
(4.11) 



~'~ 1 48 216 J TV 



(4.12) 



For the infrared factorisation of the two-loop results, the renormalised next-to-leading 
order helicity amplitude coefficients are needed through to 0(e 2 ). Their decomposition in 
colour structures is straightforward: 



nf Mite (u,v) = C F a£>(u,v). 



(4.13) 



The expansion of the coefficients through to e 2 yields HPLs and 2dHPLs up to weight 4. 
The explicit expressions are of considerable size, such that we only quote the e°-terms in 
the appendix. To this order, the coefficients had been derived previously [38, 39] in terms 
of logarithms and dilogarithms. The expressions through to 0(e 2 ) in FORM format are 
appended to the arXiv submission of this article. 
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The finite two-loop remainder is obtained by subtracting the predicted infrared struc- 
ture (expanded through to O(e )) from the renormalised helicity coefficient. We further 
decompose the finite remainder according to the colour structures as follows: 

Of Mite (u,v) = (iV 2 - 1) + ^S#> («,«)) + C F N F C%\u,v) 

+ C F N Fy D ( £\u,v) , (4.14) 

where the last term, as discussed above, is generated by graphs where the virtual gauge 
boson does not couple directly to the final-state quarks, and is different from zero only for 
the V = Z, 7* case. 

The helicity coefficients contain HPLs and 2dHPLs up to weight 4. The size of each 
helicity coefficient is comparable to the size of the helicity-averaged tree times two-loop 
matrix element for 3j production quoted in [20], and we decided not to include them 
here explicitly. The complete set of coefficients in FORM format is attached to the arXiv 
submission of this article. 



5. Checks on the result 

Several non-trivial checks were applied to validate our results: 

1. All seven tensor coefficients in (2.10) were computed. We validated that they fulfil the 
following relations, which follow from the symmetry properties of the tensor under 
an interchange p\ p 2 : 

A 2 l(si3,S 2 3,Si23) = -Ai 2 (s23,Sl3,Si23), 

-422(,Sl3,S23,Sl23) = -in(«23, $13, S123), 

-423(«13,S23,Sl23) = -Al 3 (s 2 3, «13, Sm), 

B(s 13 ,S 2 3,S 12 3) = B(s 2 3, S13, S123). (5.1) 

(2) 

2. We computed explicitly the helicity coefficients for all the diagrams in classes I§ 

(2} 

and ig , which should yield a vanishing contribution due to Furry's theorem. Each 
diagram gives a non-vanishing contribution, a full cancellation is obtained only in the 
sum of all diagrams. 

3. The IR singularity structure of our result agrees with the prediction of the Catani 
formula [37]. 

(2) 

4. We compared the helicity coefficients Q z for qq — > Z*y, with those for 7* —> qqg [21]. 
As explained in Appendix C, the unrenormalised two-loop f^- coefficients can be 
decomposed according to their colour structure as: 

ng'™> (u, v) = N*A^ (u, v) + B§>^ (u, v) + ±C§*rt («, t,) 

+ NN F D^\u,v) + ^Eg^\u,v) + N F F^\u,v) 
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+ N Fy \± - N) C$> w Hu t v) . (5.2) 

We checked that in the decay kinematical configuration Z — > qqj, prior to UV 
renormalisation and IR subtraction, the following identities are fulfilled 

B^(y,z) = - Cg*' m \ y ,z), 

C^\y,z) = -2Eg^\y lZ ), (5.3) 
D^\y lZ ) = -AG§^\y,z), 

which follow from the structure of the underlying two-loop diagrams. UV renormali- 
sation and IR subtraction of the Z — > qqj and Z — > qqg amplitudes differ. However, 
two of the above relations are unaffected by renormalisation and retain the same IR 
structure, such that we obtain: 

B^\u,v) = - Cg^ c \u,v), 

D^ c \u,v) = - AG^ jMitc \u,v) , [ ' } 

which also remain true after analytic continuation. 



6. Conclusions and Outlook 



In this paper we derived the two-loop corrections to the helicity amplitudes for the processes 
qq — > W ±r y and qq —> Z j. Our calculation was performed in dimensional regularisation 
by applying ti-dimensional projection operators to the most general tensor structure of the 
amplitude. Our results are expressed in terms of dimensionless helicity coefficients, which 
multiply the basic tree-level amplitudes, expressed in four-dimensional spinors. By applying 
Catani's infrared factorisation formula, we extract the finite parts of the helicity coefficients, 
which are independent on the precise scheme used to define the helicity amplitudes. We 
provide compact analytic expressions for the two- loop helicity coefficients in terms of HPLs 
and 2dHPLs. 

These amplitudes are relevant to the NNLO corrections to Vj production at hadron 
colliders. This process provides a direct access to the photonic couplings of the weak gauge 
bosons and is a crucial test of the structure of the electroweak theory at high energies. 
With the amplitudes derived here, in combination with the amplitudes relevant to Vjj 
production at NLO [9], all ingredients to this NNLO calculation are now available. Since 
the leading order contribution to this process does not contain any QCD partons in the final 
state, the well-established qx subtraction method [12] could be used for this calculation. 
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A. One-loop helicity amplitudes 



We list here the e -terms of the one-loop helicity coefficients for the groups of diagrams 
defined above: 



a [ a ] = C2 + G(-v, u)i 7T - G(-v, u)H(0, v) + G(-v, 0, u) + G(0, u)H(0, v) + H(0, v)i tt 
- H{0,0,v) + H(l,0,v) + O(e) , 

(l-v) 2 



G(l,u)H(l,v) 



G(l - v, u)H(0, v) + G(l, 1 - v, u) + G(l, u)i it - G(l, u)H(0, v) 
(1-u- v){2 -2v-u + 3uv- u 2 ) 



+ 



2w(l-u) 2 



% 7r + G(1 -v,u) -H(l,v) 



v(2 - 2v - 2u + 3uv) , s (7u-7-v) 
f — ~ " H(0, v) + ^— + 0(e) 



2u(l-u) 2 
(3) _ 4(1 -u-v) , 
(l-v) 



2(1 - u) 



(A.l) 



a { p = C 2 + G(-w, u)i 7T - G(-d, u)#(0, v) + G(-v, 0, u) + G(0, u)if(0, u) + (0, u)i tt 
- H(0, 0, „) + * (1, 0, „) + V( ^ 3V v) + 2 U) m v) + ( " 9 2( t^ U) + ^(e) 

G(l - u, it)i/(0, u) + G(l, 1 - v, u) + G(l,u)i tt - G(l,u)H(0, v) 

(2-2v-u + 3uv- u 2 ) 



(2) _ (l-v + uv) 
°" ~ u 2 



G(l,u)H(l,v) 



+ 



2u(l - u) 

v(2 - 4v + 2v 2 + 3uv - 3uv 2 - 2u 2 + 3u 2 v - u 3 ) 



in +G(l-v,u)-H(l,v) 



2u(l - u){\ - v) 2 



a (3) - 



4u 



(1-v) 



0(e) 



(A.2) 



,.(2) 



V(1 - U - V) H(0,v)+ { \- U -^ + O(e) 



2(1 -v) 2 

v(l — u — v) 
u 2 

G(l,u)H(l,v) 



2(1 -v) 

G(l-v,u)H(0,v) + G{l,l-v,u) + G{l,u)in - G(l,u)H(0,v) 
v(l -u-v){2-3u) 



+ 



in + G{1 - v,u) - H(l,v) 



2u(l-u) 2 

(1-u - v)(v -u) v(l — u — v)(2v - 2v 2 - 2uv + 3uv 2 - 2u 2 v + u 3 ) 



2(l-u)(l-v) 



2u{l-u) 2 {l-vf 



4 3) - 0(e) . 



H(0,v) + O(e) , 

(A.3) 



It should also be noted that these finite pieces of the one-loop coefficients can equally 
be written in terms of ordinary logarithms and dilogarithms, see [38,39]. The reason to 
express them in terms of HPLs and 2dHPLs here is their usage in the infrared counter-term 
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of the two-loop coefficients, which cannot be fully expressed in terms of logarithmic and 
poly logarithmic functions. 



B. Two-loop leading colour amplitudes 

The analytical expressions for the finite remainders of the two-loop amplitudes, as de- 
fined in (4.14), are of considerable size. For illustration, we only quote the leading colour 
contributions here. 



A™ = - 

a 4 



G(— v, 1 — v, —v,u)iir + G(—v, 1 — v, — v, u) if(0, v) — G(—v,l — v, — v, 0, u) 

- G(-v, l-v,u)( 2 - G(-v, 1 - v, u) if (0, v)iir + G(-v, l-v,u) if (0, 0, v) 

- G(-v, 1 - v,u)H(l,0,v) - G(-v, 1 - v,0,u)H(0,v) -2G(-v,-v,-v,u)iir 
+ 2 G(-v, -v, -v, u) if (0, v)-2 G(-v, -v, -v, 0, u) - 4 G(-v, -v, u) ( 2 

- 2 G(-v, -v, u) if (0, v)in + 2 G(-v, -v, u) if (0, 0, v) + 2 G(-v, -v, 0, u) i tt 

- 2 G(-v, -v, 0, u) if (0, v) + 2 G(-v, -v, 0, 0, u) - 2 G(-v, u) Q z 

- G(-v, u) i tt C 2 - 3 G(-v, u) if (0, 0, v) i tt + 3 G(-v, u) if (0, 0, 0, v) 

- G(-v, u) if (0, 1, 0, v) + G(-v, u) H(l,v) C 2 + G(-v, u) if (1, 0, v) i tt 

- G(-v, u) if (1, 0, 0, v) + G(-v, u) if (1, 1, 0, v) + 3 G(-u, 0, -v, u) i tt 

- 3 G(-v, 0, -t>, u) if (0, v) + 3 G(-v, 0, -w, 0, it) - G(-v, 0, u) ( 2 

+ 2 G(-v, 0, u) if (0, v)in-3 G(-v, 0, u) if (0, 0, u) + G(-v, 0, u) if (1, 0, v) 
+ 2 G(-u, 0, 0, u) if (0, v) + 2 G(0, 1 - v, -v, u) i it - 2 G(0, 1 - u, -v, u) if (0, u) 
+ 2 G(0, 1 - v, -v, 0, u) + 2 G(0, 1 - u, u) C2 + 2 G(0, 1 - u, u) if (0, v) i tt 

- 2 G(0, 1 - v, u) if (0, 0,v) + 2 G(0, l-v,u) H(1,0, v) + 2 G(0, 1 - v, 0, it) if (0, v) 
+ 2 G(0, -u, -v,u)in-2 G(0, -u, -u, it) if (0, v) + 2 G(0, -w, -u, 0, u) 

- 2 G(0, -v, u) C 2 + G(0, -v, u) H(0, v)iir-2 G(0, -v, u) H(0, 0, w) 
+ G(0, -u, 0, u) H(0, v) - G(0, it) H(0, v) ( 2 + 2 G(0, u) H(0, 0,v)iiT 

- 3 G(0, u) i? (0, 0, 0, v) + G(0, u) H(0, 1, 0, u) + 2 G(0, u) (1, 0, 0,v)-2 G(0, 0, -u, w) i tt 
+ 2 G(0, 0, -v, u) H(0, v)-2 G(0, 0, -v, 0, u) + 2 G(0, 0, u) H(0, 0,v)-2 H(0, v) ( 3 

- H(0, v) i tt C2 - 3 i? (0, 0, 0, v) i tt + 3 H(0, 0, 0, 0, v) - H(0, 0, 1, 0, v) + H(0, l,v)( 2 
+ H(0, l,0,v)in- H(0, 1, 0, 0, v) + H(0, 1, 1, 0, v) + 2 if (1, u) C 3 

+ 2 if(l, 0, u) C2 + 2 if(l, 0, 0, v) i tt - 2 if(l, 0,0,0,u) + 2 if(l, 0,l,0,u) + 2 if (1,1,0,0,?;) 



_5_ 11 

16 U + 12 



G(— v, —v, u) i tt + G(—v, —v, u) H(0, v) — G(—v, —v, 0, u) 



- G(-v, u) if (0, 0, v) + G(-v, 0, u) if (0, v) - G(-v, 0, 0, u) - G(0, 0, u) if (0, v) 



G(-w,«)fl-(0,t;)C2 + H(0,0,i;)C2 
v 



11 



G(-W,«)C2 



16 (u + ?;) 2 V / v ' 7 



1 — 3 1) — 3 u 

(2-2v - Au + buv + 4u 2 -5u 2 v - 2u 3 j - G(0, -u, w) ? tt 



8u(l-u) 3 
+ G(0, -u, u) if (0, ?;) - G(0, -u, 0, u) 
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+ 



(1 — u — v) 
8w(l- w) 3 
(1 — u — v) 



4u 2 



(2-2v -3u + 5uv - 5u 2 v + u 3 ^j 
1 — v + u 



-2iir( 2 - H(l,0,v)iw 



+ G(l — v,l — v, —v, u) i tt — G(l — v,l — v, —v, u) H(0, v) 

+ G(l-v,l- v, -v, 0, u) + G(l-v,l- v, u) ( 2 +G(l-v,l- v, u) H(0, v) i tt 
-G{l-v,l- v, u) H(0, 0, v) + G(l - v, 1 - v, u) H(l, 0, v) + G(l-v,l- v, 0, u) H(0, v) 
+ G(l - v, u) C 3 + G(l - w, u) H(0, v) ( 2 + G(l - v, u) H(0, 0,v)in 

- G(l - v, u) H(0, 0, 0, v) + G(l - v, u) H(0, 1, 0, w) + G(l - w, u) i? (1, 0, 0, v) 

- G(l - v,0, -v,u)iw + G(l - v,0,-v,u)H(0,v) - G(l - u,0, -u,0,u) 

+ G(l - u, 0, u) H(0, 0, w) + G(l, 1 - -u, u) iir - G(l, 1 - v, -v, u) H(0, v) 

+ G(l, 1 - v, -v, 0, u) + G(l, 1 - v, u) C2 + G(l, 1 - v, u) H(0, v) i tt 

-G(l,l-v, u) H(0, 0, v) + G(l, l-v,u) H(l, 0, v) + G(l, 1 - v, 0, u) ff(0, v) 

+ G(l,u) ( 3 -2 G(l, u) i tt C 2 - G(l, u) #(0, u) C 2 - G(l,u) H(0, 0, t>) i it 

+ G(l,u)H(0,0,0,v) - G(l,u)H(0,l,0,v) - G(l, u) H(l, v) ( 2 - G(l, u) H(l, 0, v) iir 

+ G(l, u) H(l, 0, 0, u) - G(l,u) H(l, l,0,v) + G(l, 0, -v, u)in- G(l, 0, -v, u) H(0, v) 

+ G(l, 0, -t), 0, u) - 2 G(l, 0, u) C2 - G(l, 0, u) (0, 0, v) - G(l, 0, u) H(l, 0, w) 

+ tttt tt ^ ( - v - v 2 + 2v 3 - 4uv + 5uv 2 - 4u 2 + 7u 2 v + 4u 3 ) \i7T + G(0,u) 

16 (1 — u) (u + vr V / L 



+ 



16 (1 - u) {u + v) 
1 

144 (1 - u) 2 (1 - u - v) 

G(0, u) H(0, v) + H(0, v)iir 



313 - 331 v + 18 v 2 - 975 u + 716 u v - 54 u v 2 + 1011 u 2 



385 u 2 v - 349 u 3 



+ 



36 - 108 v + 108 v 2 - 36 « 3 + 169 u + 29 u v ~ 270 u 



H ( -8v + 8v 2 -9u + 25uv+l7u 2 ) 

16 u (u + v) V / 

+ 1 f 

144u(1-m) 2 (1-m-w) V 

+ 72 u v 3 - 759 u 2 + 320 it 2 v + 126 u 2 w 2 + 867 u 3 - 241 u 3 v - 313 u 4 ) i? (1, 0, u) 

(72 - 216 v + 216 v 2 - 72 v 3 + 25 u + 389 u v - 558 u w 2 

144 u(l - u) 2 (1 - u - v) V 

144 at 3 - 543 u 2 - 76 u 2 w + 306 u 2 w 2 + 723 u 3 - 97 u 3 v - 277 u 4 ) ( 2 

— f-36i; + 72w 2 -36w 3 -313u + 475uw- 234ui; 2 



+ 



144 u (1 - u) 2 (1 - it - t>) 
+ 72 u w 3 + 975 u 2 - 896 u 2 v + 198 u 2 v 2 - 1011 u 3 + 457 u 3 v + 349 u 4 ) - G(-u, u) i tt 

+ G(-w, u) H(0, v) - G(-v, 0, u) + i? (0, 0, u) 

+ 7, — 77^ nr ( 4 - 6 v + 2 v 2 - 10 u + 16 u v - 5 u w 2 + 6 u 2 - 16 u 2 v + 5 u 2 v 2 

8u (1 — m) 3 V 

+ 2 u 3 + 6 u 3 1) - 2 u 4 ) G(l - w, -u, u) i tt - G(l - v, -v, u) H(0, v) + G(l - v, -v, 0, u) 
+ G(l - v, u) ( 2 + G(l - v, u) H(0, v) i tt - G(l - u, it) i?(0, 0, u) + G(l - v, u) H(l, 0, u) 
+ G(l-u,0,u)fl'(0,u) 

- - / f6 - 12 v + 6v 2 - 4?i + 30uw - 15uv 2 - 24 u 2 - 30u 2 v + 15u 2 v 2 

12 u (1 - u) A \ 
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+ 36u 3 + 12u 3 v - 14 u 4 ) G(0,u)C 2 

+ 7T, — r- ^ f6 - 12 1; + 6« 2 + 7?i + 30uw - Ihuv 2 - 57 u 2 - 30u 2 v + 15u 2 v 2 

24 u (1 - u) 3 \ 

+ 69 u 3 + 12 u 3 v - 25 u 4 ) - G(0, u) if (1, 0, v) - if (1, v) ( 2 - if (1, 1, 0, v) 

H J- ( 12 - 18 v + 6 v 2 - 8 u + 48 u v - 15 u v 2 - 48 u 2 - 48 u 2 v + 15 u 2 V 

24 u (1 - u) 3 V 

+ 72u 3 + 18u 3 w-28u 4 )c3 



+ 



12-l8v + 6v 2 + Uu + 48uv-15uv 2 -lUu 2 -48u 2 v + I5u 2 v 2 



24u(l-u) 3 
+ 138 u 3 + 18 it 3 u - 50 u 4 , ) if (1, 0, 0, u) 

+ — T- ^ f - 6w + 6w 2 - 44m + 12 mi; - 15?iv 2 + 132 u 2 - 12m 2 w + 15u 2 t> 

24 u (1 - u) 3 \ 



132u 3 + 6u 3 i; + 44u 4 



+ 24m (1 - u) 
- 66 u 3 + 6 u 3 v + 22 u 4 



1 22 u -\- 12 uv — 15 uv 2 -\- 66 u 2 — 12 u 2 v 15 u 2 u 2 



- G(0, it) if (0, 0, v) + if (0, 0, 0, v) 

huV 3 (- 6v+6v 

-H(0,0,v)in- H(0,l,0,v 

+ 7T< — T- f 6 v - 6 v 2 - 22 u - 12 u v + 15 u v 2 + 66 it 2 + 12 u 2 v - 15 u 2 v 

24 u (1 - u) 3 V 



i (1 - u) 3 
- 66 u 3 - 6 u 3 v + 22 u 



)h(0,v)( 2 , 



.C9i 11 . . 49 , 7 



+ 



(1-^) S 
4 m 2 



13 35 22 

— G(l, u) if (0, v)in-- G(l, u) if (0, l,v) + — G(l - «, u) if (0, 0, v) 
6 6 3 
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- — G(l, u) if (0, 0,v) + 2 G(l -t),l-t),«)C2 + 2G(l-t),l- v, u) if (0, v) i it 

-2G(l-v,l-v,u)H(Q,l,v) -2G(l-v,l - v, 1, 1 - v, u) - 2 G(l - v, 1 - v, 1, u) iir 

+ 2G(l-v,l-v,l,u) if (0, v) + 2G(l-v,l- v, 1, u) if (1, u) + 2 G(l - u, 0, 1 - v, u) if (0, v) 

+ 2 G(l - t), 0, 1, 1 - d, u) + 2 G(l - u, 0, 1, u) i tt - 2 G(l - v, 0, 1, u) if (0, v) 

- 2 G(l - w, 0, 1, u) if (1, u) + G(l - v, 1, 1 - u) (0, v)-A G(l - u, 1, u) Q 2 

- G(l - v, l,u)H(0,v)iir+ G(l - w, l,u) H(0, l,v) - G(l - v, l,u) H(l,0,v) 
+ 2 G(l - v, 1, 1, 1 - t>, u) + 2 G(l - u, 1, 1, u) i it - 2 G(l - v, 1, 1, u) H(0, v) 

- 2 G(l - u, 1, 1, u) H{1, v) + 2 G(l, 1 - v, 1 - u, u) H(0, v) - 4 G(l, 1 - u, u) C2 

- G(l, 1 - u, u) if (0, u) i 7T + G(l, l-v,u) if (0, 1, u) - G(l, 1 - u, u) if (1, 0, u) 
+ 3 G(l, 1 - u, 1, 1 - v, u) + 3 G(l, 1 - v, 1, u) i tt - 3 G(l, 1 - v, 1, u) ff(0, u) 

- 3 G(l, 1 - v, 1, u) ff(l, w) - 3 G(l, u) Ca - G(l, u) i tt C2 + 4 G(l, u) ff(l, v) (2 
+ G(l, u) if (1, 0,v)in- G(l, u) if (1, 0, 1, v) + G(l, u) if (1, 1, 0, u) 

- G(1,0,1- u,u)ff(0,u) - G(l,0,l,l-u,u) - G(l,0,l,u)iir + G(l,0,l,u)H(0,v) 
+ G(l,0,l,u)ff(l,w) + 2G(1,1,1- v,l- v,u) + 2G(l,l,l- v,u)iTr 

- 2G(1, 1, 1 - v,u)H(0,v) - 2G(1, 1, 1 - u, u)if(l,w) - 2G(1, (2 
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+ 



2 G(l, 1, u) H(l, v) i tt + 2 G(l, 1, u) H(l, 0, v) + 2 G(l, 1, u) H(l, 1, u) 

2 G(l, 1, 1, 1 - u, u) - 2 G(l, 1, 1, u) i tt + 2 G(l, 1, 1, u) H(0, v) + 2 G(l, 1, 1, u) H{1, v) 

(2-2tJ-2w + 3iiuj i?(0, 0, v) 



llv 



12u(l - ?i) 2 
(1 — u — v 



+ 7^r — % (2-2v-u + 3uv-u 2 ) \ -G(l-v,u)H(l,0,v) + H(l,0,v)iiT 
8u (1 — m) 2 V / L 

-if(l,0,M) + #(1,1,0,*;) 

HzJiZ ^ f6-6i;-3u + 5uw-3u 2N ) [g(1 - t>, u) i tt - H(l, v) i tt 
16 m (1 — u) z \ / L 

1 (22 - 22w-llu + 39?iv-llw 2 ) J ff(l,0,u) 



+ 



+ 



24u(l-u) 2 
(1 — u — v) 



26 - 26 v - 13 u + 51 u v - 13 u 2 



-z^ u (1 — u) 2 . 
+ G(l-u,u)i2"(l,u) -H(l,l,v) 
(1 — u — v) 



— G(l — v, 1 — u,u) 



+ 



( - 554 + 288 ( 2 + 554 v - 288 v ( 2 + 433 u - 144 u C 2 - 909 u v 



288u(l-u) 2 

+ 432u«C 2 + 121 u 2 - 144w 2 C 2 ) H(l,v) 

(1 — u — v) 
288 u(l- u) 2 
(1 - u - v) / 
8u 2 (l -u) 2 V 
1 



(554- 554w-433u + 909uw- 121 u 2 ) G(l-v,u) 



3 — 3v — u + 2uv — u +3u w — u" 



G{l,l-v,u)in- G(l,u) H(l, v) i tt 



(5184 - 5184 v - 87239 u - 4788 uv + 82055 u 5 



10368 u(l-u) 

H — ^ — (2-8v + 6v 2 -3u+l0uv-9uv 2 -2u 2 v + u 3 ) [g(1 - u, 1, 1 - u, u) 

8u (1 — w) 2 V / L 

+ G(l - u, 1, u) i tt - G(l - v, 1, u) ff(0, t>) - G(l - v, 1, u) ff(l, t>) 

+ — t — (2-2v 2 -3u + 2uv + 3uv 2 -2u 2 v + u 3 ) [g(0, 1 - w, u) #(0, u) 

8u (1 — m) V / L 

+ G(0, l,l-v,u) + G(0, 1, u) i tt - G(0, 1, u) H(0, v) - G(0, 1, u) u) 

- (aA- 114i; + 70w 2 -66m+158uw- 105uv 2 - Uu 2 v + 22u 3 ) H(0,1,v) 

48 m (1 — u) z \ J 

+ T7, — r- 77 ( 44 - 70 v + 26 v 2 - 66 u+ 114 uu - 39uv 2 - AAu 2 v + 22u 3 ) H(0,v)in 

48 u (1 - u) 2 V / 

- 777-. — 7z 77 (480 C 2 + 108 C 3 - 882 v( 2 - 216 v ( 3 + 402 v 2 ( 2 + 108 v 2 ( 3 - 660 u ( 2 

144 u (1 — uy V 

- 551 u C 3 + 1326 u u C2 + 324u v C3 - 567 u!) 2 ( 2 - 162 u v 2 ( 3 - 120 u 2 C 2 + 778 u 2 ( 3 

- AAA u 2 v ( 2 - 108 u 2 v ( 3 + 300 u 3 ( 2 - 335 u 3 Q 3 ) 

^ ^ f - 1662 + 216 C 2 + 3324 w - 432 v ( 2 - 1662 w 2 + 216 v 2 ( 2 + 5324 u 

864 u (1 — u) z \ 

+ 270uC 2 - 5292 u v + 648 u v ( 2 + 2727uw 2 - 324uw 2 C 2 - 5662 u 2 - 1188 u 2 C 2 + 1968 u 2 v 

- 216 u 2 v ( 2 + 2000 u 3 + 702 u 3 ( 2 ^ i tt 
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864 u(l - u) 
+ 2363 m 3 ) H(0,v 



- f - 1662 v + 1662 v 2 + 2363 u + 2526 uv- 2727 uv 2 - 4726 u 2 - 864 u 2 w 

1 — uY \ 



+ o o / 1 ^ f3-6M + 3i; 2 -6u+14Mi;-8uw 2 + 3u 2 -8u 2 w + 6?i 2 

1 — u r v 



8ii 2 (l-u) 2 
x G(l -M,u)H(0,w)i7r 



» 2 ) 



- — (^13 — 26 v -\- 13 v 2 — 32 u -\- 70 uv — 38 u v 2 -\- 22 u 2 — 50 u 2 v -\- 31 u 2 v 2 

1 — uV V 



24u 2 (l -u) 



+ 

+ 6« 3 ^3w 4 ) -G(l-v,l-v,u) H(0, v) - G(l, 1 - v, 1 - v, u) + G(l, 1 - v, u) H(l, v) 
-G(l,u)H(l,l,v) 

+ ^ o } ^ f 13 - 26 w + 13 v 2 - 26 u + 46 u w - 20 u v 2 + 13 u 2 - 20 u 2 v + 4 u 2 w 2 ) 

24 u 2 (1 - u) 2 V / 

x G{l-v,u)H(0,l,v) 
1 



+ 



13 - 26 u + 13 w 2 - 20u + 34uv- Uuv 2 + 4u 2 - 2u 2 v - 5u 2 v 2 - 6u 3 v 



24u 2 (l-u) 2 

+ 3« 4 ) - G(l, 1, 1 - v, u) - G(l, l,u)in+ G(l, 1, u) ff(0, u) + G(l, 1, u) u) 

- o / ^7 fl3 - 26 w + 13 w 2 - 8u+ 10 um - 2uw 2 - 14 u 2 + 34u 2 t; - 23u 2 v 2 

24 u 2 (1 - u) 2 V 

-18u 3 v + 9u 4 ) G(1-w,m)C 2 



+ 



22 — 44 v -\- 22 v 2 — 50 u -\- 106 uv — 56 u v 2 -\- 31 u 2 — 68 u 2 v -\- 40 u 2 v 2 



24u 2 (l-u) 

+ 6« 3 «-3!i 4 j G(l, 1 - v, u) H(0, v) - G(l,u) H(l, 0, v) 

67 — 134 v -f- 67 v 2 — 128 u -\- 250 uv — 122 u v 2 -\- 58 u 2 — 110 u 2 v -\- 49 u 2 v 2 



24u 2 (l-u) 2 
-6ti 3 u + 3u 4 j G(1,m)C 2 

+ o\ ^ f277 — 554 v -\- 277 w 2 - 650 u + 1264 o- 614 u v 2 + 430 u 2 - 770 u 2 v 

144 u 2 (1 - u) 2 V 

f G(l, 1 - v, u) + G(l, u) i 7T - G(l, it) ff(0, v) 



+ 340 u 2 u 2 - 18 m 3 + 60 u 3 v - 39 w 4 



-G(l,u)H(l,v) 



+ 



1 (■ 

fl-w) 2 V 



277 — 554 v -\- 277 v 2 - 518 u + 922 u v - 404 u w 2 + 232 u 2 - 296 u 2 v 



144u 2 (l -u) 

+ 25u 2 w 2 - 18m 3 - 72m 3 m + 27m 4 ) G(l - v, u) H(0, v) , 



A (3) - 



(1 — u — v) 
2759 . 2759 



5 11 11 49 7 7 389 

^C 2 - T 7,-C 2 - T 7,^(0,,)C 2 + ^C4 + i ^C3 + ^(0,,)C3 + T ^C3 



864 



■l 7T ■ 



864 



ir(o,«)- 



81659 
10368 



a< 4 >=o. 



(B.l) 
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= i — G(— v, 1 — v, —v,u)itt + G(—v, 1 — v, — v, u) H (0, v) — G(—v,l — v, — v, 0, u) 

- G(-v, l-v,u)( 2 - G(-v, 1-v, u) H(Q, v)in + G(-v, l-v,u) H(0, 0, v) 

- G(-v, 1 - v, u) H(l, 0, v) - G(-v, l-v,0, u) H(0, v)-2 G(-v, -v, -v, u) i tt 
+ 2 G(-v, -v, -v, u) H(0, v)-2 G(-v, -v, -v, 0, u) - 4 G(-v, -v, u) ( 2 

- 2 G(-v, -v, u) H(0, v)in + 2 G(-v, -v, u) H(0, 0, v) + 2 G(-v, -v, 0,u)iir 

- 2 G(-v, -v, 0, u) H(0, v) + 2 G(-v, -v, 0, 0, u) - 2 G(-v, u) (3 - G(-v, u)iir( 2 

- 3 G(-v, u) H(0, 0, v) i tt + 3 G(-v, u) H(0, 0, 0, v) - G(-v, u) H(0, 1, 0, v) + G(-v, u) H(l, v) Q 2 
+ G(-v, u) H(l, 0,v)iir- G(-v, u) H(l, 0, 0, v) + G(-v, u) H(l, 1, 0, v) 

+ 3 G(-v, 0, -v,u)iir-3 G(-v, 0, -v, u) H(0, v) + 3 G(-v, 0, -v, 0, u) - G(-v, 0, u) C2 
+ 2 G(-t>, 0, it) ff(0, v) i tt - 3 G(-t>, 0, tt) ff(0, 0, v) + G(-t>, 0, u) H (1, 0, u) 
+ 2 G(-d, 0, 0, u) i? (0, u) + 2 G(0, 1 - -v, tt) i tt - 2 G(0, 1 - v, -v, u) H(0, v) 
+ 2 G(0, 1 - v, -v, 0, tt) + 2 G(0, 1 - u, tt) C2 + 2 G(0, 1 - u, tt) H(0, v) i tt 

- 2 G(0, 1 - v, tt) H(0, 0,v) + 2 G(0, l-v,u) H(1,0, v) + 2 G(0, 1 - v, 0, tt) ff(0, v) 

+ 2 G(0, -u, -v,u)i-K-2 G(0, -t>, -u, it) H(0, v) + 2 G(0, -w, -u, 0, u) - 2 G(0, -v, u) ( 2 
+ G(0, -v, u) H(0, v)iir-2 G(0, -v, u) H(0, 0, w) + G(0, -u, 0, u) H(0, v) - G(0, u) H(0, v) ( 2 
+ 2 G(0, u) i? (0, 0, w) i it - 3 G(0, u) H(0, 0, 0, u) + G(0, u) ff(0, 1, 0, v) 
+ 2 G(0, u) H(l, 0, 0, u) - 2 G(0, 0, -v, u) i tt + 2 G(0, 0, -w, u) H(0, v) 

- 2 G(0, 0, -w, 0, u) + 2 G(0, 0, u) H(0, 0, u) - i? (0, v) i tt ( 2 

- 3 H(0, 0, 0, u) i tt + 3 ff(0, 0, 0, 0, u) - H(0, 0, 1, 0, u) + H(0, 1, u) C2 

+ #(0, 1, 0, v) i tt - H(0, 1, 0, 0, u) + /f(0, 1, 1, 0, v) + 2 H(l, v)( 3 + 2 H{1,0, v) ( 2 

+ 2 H (1, 0, 0, v) i tt - 2 H(l, 0, 0, 0, v) + 2 0, 1, 0, v) + 2 H(l, 1, 0, 0, u) 
11 



+ 



12 



G(-w, -u, m) i tt + G(-v, -v, u) H(0, v) - G(-v, -v, 0, u) - G(-v, u) H(0, 0, v) 



+ G(-v, 0, u) H(0, v) - G(-v, 0, 0, u) - G(0, 0, tt) H(0, v) 



39 _ 5 
+ 32 C4+ 4 

1 

+ 



G(-v, u) H(0, v) ( 2 + H(0, v) C 3 - H(0, 0, v) C2] + \i n Cs - y u) C2 



8u(l-u) 2 (l-f) 



— f 2 - 4 1) + 2 w 2 - 8 u + 13 tt v - 5 u v 2 + 12 u 2 - 16 u 2 v + 5 m 2 w 2 
vV V 



- 5 u 3 + 3 u 3 1) + tt 4 ) - G(0, -t), u)iTT + G(0, -u, tt) ff(0, v) - G(0, -u, 0, u) 

+ ^ 4 ^ ^1 - v + u^ G(l - u, 1 - t;,-t;,u)t7r- G(l - v, 1 - v, -u, tt) H(0, v) 

+ G(l-v,l- v, -v, 0, u) + G(l-v,l- v, u) ( 2 +G(l-v,l- v, u) H(0, v) i tt 
-G(l-v,l- v, u) H(0, 0, v) + G(l-v,l- v, tt) H(1,0, v) +G(l-v,l- v, 0, tt) H(0, v) 
+ G(l - v, tt) Cs + G(l - v, tt) H(0, v) ( 2 + G(l - u, tt) H(0, 0,v)in- G(l - t>, tt) if(0, 0, 0, u) 
+ G(l - v, tt) H(0, 1, 0, t>) + G(l - u, tt) i?(l, 0, 0, t>) - G(l - v, 0, -t>, tt) i tt 
+ G(l - u, 0, -u, tt) ff(0, t>) - G(l - v, 0, -v, 0, tt) + G(l - v, 0, tt) i? (0, 0, v) 
+ G(l,l-v, -v, tt) i tt - G(l, 1 - u, -t>, tt) i? (0, u) + G(l, 1 - t>, -u, 0, m) + G(l, 1 - v, tt) C2 
+ G(l,l-v, tt) ff(0, d) i tt - G(l, 1 - t>, tt) if(0, 0, u) + G(l, 1 - v, tt) 0, v) 
+ G(l, 1 - v, 0, tt) #(0, v) + G(l, u) Cs - 2 G(l, «) i tt ( 2 - G(l,u) H(0, v) ( 2 
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- G(l,u) H(0, 0,v)iir + G(l,u) H(0, 0, 0, v) - G(l, u) H(0, 1, 0, v) - G(l, u) H(l, v) ( 2 

- G(l,u) H(l,0,v)in + G(l,u) H(l,0,0,v) - G(l,u) H(l, 1,0, v) + G(1,0, -v,u) in 

- G(l, 0, -v, u) H(0, v) + G(l, 0, -v, 0, u) - 2 G(l, 0, u) Q 2 - G(l, 0, u) H(0, 0, v) 

- G(l,0,u)H(l,0,v) 

- 7777 T7 77 (3v- 19w 2 + 16w 3 - 17 uv + 36uv 2 - u 2 + 39u 2 v + 19u 3 ) G(0,u) 

48 (1 - v) (u + v) 2 \ ) v ' 

+ ttttt; \-, 77 ( - 54 v - 2813 v 2 + 2867 v 3 - 6004 uv + 6058 uv 2 - 3137 u 2 

864(1 - v)(u + v) 2 V 

+ 3245 u 2 w + 54 u 3 ) z tt 

+ ttttt; 777 77 f 54 v - 3533 v 2 + 4537 v 3 - 1058 w 4 - 6580 uv + 8822 u v 2 

864 (1 — (u + v) z \ 

- 1753 uv 3 - 3155 u 2 + 4735 u 2 v - 602 u 2 v 2 + 396 u 3 + 93 u 3 v ) H(0, v) 

+ — — f 295 - 912 v + 939 v 2 - 322 v 3 - 626 u + 1700 u v 

144 (1 - u) (1 - v) 2 (1 - u - v) V 

- 1432 Mt) 2 | 358 u v 3 + 313 u 2 - 761 u 2 v + 421 u 2 v 2 + 18 u 3 + 9 u 3 v ) #(0, u) i tt 



( 



295 - 714 1; + 543 v 2 - 124 v 3 - 626 u + 1238 u w 



144 (1 - u) (1 - -u) 2 (1 - w - v) 

- 772 u v 2 + 160 uv 3 + 313 ?i 2 - 431 u 2 v + 157 u 2 v 2 + 18 u 3 - 57 u 3 v ) G(0, u) £f(0, u) 

+ „o„o 7; 1 w r f - 5184 v + 5184 w 2 - 5832 u - 79967 «»+ 85799 u v 2 

10368 ?i (1 — v) (u + v) \ 

- 85799 u 2 + 95771 u 2 v + 9972 u 3 ) 

H — - — — f 36 v - 72 v 2 + 36 v 3 + 295 u - 644 u v + 421 u v 2 

144 u (1 - u) (1 - v) (1 - u - v) V 

- 72 u v 3 - 626 u 2 + 957 u 2 v - 349 u 2 w 2 + 313 u 3 - 313 u 3 v + 18 u 4 ) G(-v, u) i tt 

- G(-v, u) H(0, v) + G(-v, 0, u) 

+ — - ^ , (36- 144w + 216w 2 - 144w 3 + 36w 4 + 187u-318uw 

144 u (1 - u) (1 - v) 2 (l-u-v) V 

+ 3uw 2 + 200ui; 3 -72uw 4 -518u 2 + 914M 2 w-448u 2 w 2 + 52u 2 i; 3 + 277M 3 - 359u 3 i; 



+ 121 u 3 v 2 + 18m 4 - 57u 4 w) H(l,0,v) 



+ , , — - NO f 72 - 288 v + 432 v 2 - 288 v 3 + 72 v 4 + 19 u + 258 u v 



144 u (1 - u) (1 - v) 2 (l-u-v) 

- 717 at) 2 ! 584 uv 3 - 144 ti» 4 - 290 u 2 + 290 u 2 v + 116 u 2 v 2 - 116 it 2 w 3 + 181 u 3 

- 167 u 3 v + 25 u 3 v 2 + 18 u 4 - 57 u 4 v ) C2 

+ — - N „ , f - 36 « + 108 w 2 - 108 « 3 + 36w 4 - 295u + 543uw 

144 u (1 - u) (1 - v) 2 (l-u-v) V 

- 273uw 2 + 97uw 3 - 72uw 4 + 626u 2 - 659 u 2 w - Uu 2 v 2 + A7 u 2 v 3 - 313 it 3 -34u 3 w 



+ 215 u 3 v 2 - 18 u 4 + 150 u 4 v J i? (0, 0, v) 

^ -tt (2-2v-3u + 5uv-5u 2 v + u 3 ) H(l, 0, w) i tt 

1 — ii ) V / 



8u(l - u) 



- 23 - 



J- -r (6-6v + 2u + 15uv-22u 2 -15u 2 v + Uu 3 ) G(0,u)( 2 

12 u (i - uy \ j 

-t — ( 8 - 8 v - u + 20 u v - 22 u 2 - 20 u 2 v + 15 u 3 ) i n ( 2 

16 u (1 — u) A \ ) 

+ J- - (6-6v + I3u + 15uv-Uu 2 - 15u 2 v + 25u 3 ) \ - G(0, it) H(l, 0, v) 

24 u (1 - u) z \ / L 

-H(l,v)( 2 -H(l,l,0,v) 



+ 



— ^— — f 4 - 10 v + 8 v 2 - 2 v 3 - 6 u + 14 u u - 13 u v 2 + 5 u v 3 + 2 u 2 w 

tt (1 — w) 2 (1 — v) 2 V 

+ 4 u 2 u 2 - 5 u 2 v 3 + 2 u 3 - 9 u 3 v + 5 u 3 v 2 + u 4 v ) G(l - v, -v, u) i n 

- G(l - v, -v, u) H(0, v) + G(l - v, -v, 0, it) + G(l - v, u) ( 2 + G(l - v, u) H(0, v) i n 

- G(l - v, it) H(0, 0, v) + G(l - v, it) H(1,0, v) + G(l - v, 0, u) H(0, v) 
1 



+ 



(l2 - 30 v + 24 v 2 - 6 v 3 + 26 u - 46 u i> + 5 it v 2 + 15 u v 3 - 88 u 2 



24tt(l - u) 2 (1 - u) 2 
+ 182 u 2 v - 76 u 2 v 2 - I5u 2 v 3 + 50 it 3 - 115 u 3 v + 59 u 3 v 2 + 3 it 4 v ) ff(l, 0, 0, v) 

H ^ -tt (6t> - 12 w 2 + 6v 3 + 22 u- 68ui> + 61 uv 2 - 15uv 3 - Uu 2 

24 m (1 - u) 2 (1 - v) 2 V 

+ 124 u 2 v - 92 u 2 t; 2 + 15 u 2 v 3 + 22 u 3 - 59 u 3 v + 31 u 3 v 2 + 3 u 4 v ) \h(0, 0,v)in 



+ 



H(0,l,0,v) 
1 



(& v - 12 v 2 + 6 v 3 + 44 u - 112 u v + 83 u v 2 - 15 u v 3 - 



88 u 2 



24m (1 - u) 2 (1 - w) 2 

+ 212 u 2 v - 136 u 2 v 2 + 15 u 2 v 3 + 44 u 3 - 103 u 3 v + 53 u 3 v 2 + 3 u 4 v ^ G(0, u) if (0, 0, t>) 

- ff(0,0,0,u) 

H ^ -tt (l2v- 2Av 2 + I2v 3 - 77u+l06uv + uv 2 - 30uv 3 + 154 u 2 

48 u (1 - u) 2 (1 - v) 2 \ 

- 236 u 2 v + H8 u 2 v 2 + 30 u 2 v 3 -77u 3 + 124 u 3 v - 59 u 3 v 2 + 6 it 4 v ) if (0, w) C2 



+ 



-m) 2 (1-i;) 2 V 



72 - 180 v + 144 v 2 - 36 v 3 + 413 u - 790 u v + 287 u v 2 + 90 it v 3 



U4u(l-u) 2 {l-v) 

- 1042 u 2 + 2120 u 2 v - 970 u 2 v 2 - 90 u 2 v 3 + 557 u 3 - 1204 u 3 v + 611 u 3 v 2 + 18 u 4 v ) Cs, 



Al 2) = 



+ 



11 u 



12u(1-m)(1-w) 



— ^ ^2 - 4w + 2v 2 + 3uv - 'iuv 2 - 2u 2 + 3u 2 v - u 3 ^j H(0,0,v) 



1 ( 

(l-v)\ 



288u(l-v) 



144- 144w+ 133m- 133 u v + 277 



+ - — 7^ r (2-2v-u + 3uv-u 2 ) \- 3G -in (2 - G(l - v,u)H(l,0,v) 

8u (1 — u) V / L 

+ 4H(l 7 v) (2 + H(l, 0,v) in — H(l, 0, 1, v) + H(l, 1,0, v) 

H f 6 — 6w — 3u + 5uv — 3u 2 \ G(l — v,u) in — H(l,v) in 

lou (1 — tt) V / L 



- 24 - 



1 - u) ( 



24u(l - u) 



22- 22v - 11 u + 39uv - 11 u 



)H(l,0,v) 



+ 



— - (26- 26 v- 13u + 5luv- 13 u 2 ) \- G(l - v,l -v,u) 

\u (1 — u) V / L 



+ G(l-u,u)i?(l,t;) -H(l,l,v) 
1 



+ 



— (554- 1108i; + 554 v 2 - 397u + 1306m v- 909uv 2 - 271 it 2 
in V 



288«(l-u)(l- v) 
+ 157u 2 -u + 114u 3 ) G(l-v,u)-H(l,v) 



+ 



— (554- 
v) V 



1 108 v + 554 v 2 - 265 u + 1 174 u v - 909 u v 2 - 271 w 2 



288u(l-u)(l- w) 
+ 289 ?i 2 v - 18 u 3 ) i 7T 

H 7 T7 To f 2 - 10 1) + 14 v 2 - 6 w 3 - u + 8 u v - 16 u w 2 + 9 u v 3 - u 2 + 4 u 2 v 

8u (1 — u) (1 — w) 2 V 

-4« 2 o 2 + !i 3 i)j G(l - u, 1, 1 - v, u) + G(l - v, 1, u) i 7T - G(l - v, 1, u) #(0, u) 
-G(1-u,1,m)^(1,u) 



+ 



t )fl-w) 2 ( 



- u) (1 - v) 



u 3 v 



2-2v -2v* + 2v 6 -u + 2uv + 2uv z -3uv 6 -u-" + 2u z v 



2 „.2 



) G(0, 1 - w, u) i? (0, u) + G(0, l,l-v,u) + G(0, 1, u) i?r - G(0, 1, u) H(0, v) 
-G(0,l,u)H(l,v) 

7 Ki ^ f 44 - 158 w + 184 w 2 - 70 « 3 - 22 u+ 110 it v- 193 uv 2 + 105 m« 3 

48 u (1 - u) (1 - w) 2 V 

- 22 u 2 + 70 u 2 v - 61 u 2 w 2 + 13 it 3 v ) ff(0, 1, u) 

H ; ^— -r (44- 114i; + 96w 2 - 26 v 3 - 22u+ 110 - 127 uv 2 + 39uw 3 

48 u (1 - u) (1 - v) 2 \ 

- 22 u 2 + 26 u 2 v + 5 u 2 v 2 - 9 u 3 v ) £f(0, u) i tt 



— f 

- u) (1 - v) 2 \ 



160 - 454 v + 428 v 2 - 134 w 3 - 80 u + 388 u v - 497 u v 2 



48u(l-?i)(l--y) 2 
+ 189uw 3 - 80u 2 + 134 u 2 w - 41u 2 w 2 - 13u 3 v) C2 

7 ^"TT ^ ( - 554 w+ 1108-y 2 - 554 v 3 - 132 u + 384 u v - 1161 u v 2 

288u(l-u)(l- v) 2 V 

+ 909 u v 3 + 494 u 2 v - 657 ?i 2 -y 2 + 132 it 3 + 31 u 3 v ) ff(0, v) 

+ ^2 ( X - v + u v ) [ y G ( X ' u ) ^(°' w) i tt - y G(l, u) ff(0, 1, «) + y G(l - «, u) fr(0, 0, u) 
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G(l, u) H(0, 0,v) + 2 G(l - w, 1 - u, u) C2 + 2 G(l - v, 1 - u, w) H(0, v) i tt 



-2G(1-!),1- v,u) H(0, l,v) -2G(l-«,l-u,l,l- 

- 2G(1 - v, 1 - v, l,u) i-K + 2G(1 - v, 1 - v, l,u) H(0, v) + 2G(l-v,l- v, l,u) H(l, v) 
+ 2 G(l - v, 0, 1 - v, u) i? (0, u) + 2 G(l - u, 0, 1, 1 - v, u) + 2 G(l - v, 0, 1, u) i tt 

- 2 G(l - v, 0, 1, u) H(0, v)-2 G(l - w, 0, 1, u) i? (1, v) + G(l - u, 1, 1 - v, u) H(0, v) 

- 4 G(l - v, 1, u) C2 - G(l - w, 1, u) ff(0, w)i7r + G(l — v, l,u) H(0, l,v) 
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- G(l -v,l,u) H(l, 0, v) + 2 G(l - v, 1, 1, 1 - u, tt) + 2 G(l - v, 1, 1, tt) i tt 

- 2 G(l - t), 1, 1, u) H(0, v)-2 G(l - v, 1, 1, it) (1, u) + 2 G(l, 1 - t>, 1 - v, u) H(0, v) 

- 4 G(l, 1 - u, tt) C 2 - G(l, 1 - v, tt) #(0, u) i tt + G(l, 1 - v, tt) H(0, 1, u) 
-G(l,l-w,u)i?(l,0,w) + 3G(l,l-u,l,l-u,u) + 3G(l,l-t;,l,u)i7r 

- 3 G(l, 1 - «, 1, tt) if (0, v)-S G(l, 1 - v, 1, tt) if (1, v) - 3 G(l, tt) Ca - G(l,u)iirCa 
+ 4 G(l, u) if(l, v) ( 2 + G(l,u) if(l, 0,v)iir- G(l,u) H(l, 0, 1, v) 

+ G(l,u) H(l, 1,0, v) - G(l, 0, 1 - i), tt) if(0, v) - G(l, 0, 1, 1 - u, tt) - G(l, 0, l,u)in 
+ G(l, 0, 1, tt) H(0, v) + G(l, 0, 1, tt) if(l, v) + 2 G(l, 1, 1 - u, 1 - v,u) 

+ 2 G(l, 1, 1 - v, u) i tt - 2 G(l, 1, 1 - v, u) H (0, v) - 2 G(l, 1, 1 - v, u) H (1, v) - 2 G(l, 1, tt) ( 2 

- 2 G(l, 1, u) H(l, v) i tt + 2 G(l,l,tt)if(l, 0,v) + 2 G(l,l,u) H(l,l,v) 

-2G(l,l,l,l-v,u)-2G(l,l,l,u)i7T + 2G(l,l,l,u)H(0,v) + 2G{l,l,l,u)H(l,v) 

H — ( 3 - 3 v - u + 2 u v - u 2 + 3 u 2 v - u 3 ) |g(1, 1 - v, tt) i tt - G(l, tt) H(l,v) i tt 

8tr (1 — u) \ / L 

H — ^ f 13 - 13 v - 19 tt + 38 tt v + 3 tt 2 - 31 u 2 v + 3 u 3 ) 

24 u 2 (1 - u) \ J L 

- G(l - u, 1 - t>, tt) if(0, u) - G(l, 1 - t;, 1 - v, tt) + G(l, 1 - v, tt) if(l, v) - G(l, u) H(l, 1, u) 
1 



+ 



24tt 2 (l-tt) 



13-13t;-7tt+14ttt;-3u 2 + 5tt 2 t;-3u 3 ) - G(l, 1, 1 - v, tt) 



- G(l, l,u)iTT + G(l, l,u) #(0,v) + G(l, l,tt) fr(l,«) 

+ =4 7 f 22 - 22 1; - 28 u + 56 uv + 3u 2 - 40u 2 v + 3u 3 ) [g(1, 1-v,u) H(0,v) 

24 tr (1 — u) \ / L 

- G{l,u)H(l,0,v) 

J- r f67- 67t;-61u+122ut;-3tt 2 -49u 2 t;-3u 3N ) G(l,tt)C 2 

24 tt 2 (1 — u) \ J 



+ 



1 



277 - 554 v + 277 1; 2 - 355 u + 969 tt v - 614 tt f 2 + 21 tx 2 



144tt 2 (l-tt)(l-t;) 
379 u 2 v + 340 tt 2 -y 2 + 75 tt 3 - 39 tt 3 v - 18 tt 4 ) G(l, 1 - v, u) + G(l, tt) i tt 



- G(l,u) H(0, v) - G(l, tt) H(l,v) 
+ 



—z-. l —, -tt f3-9w + 9w 2 -3t; 3 -3ti+14ttt;-19ttt; 2 + 8tit; 3 -3ti 2 t; 

tt^ (1 — tt) (1 — vy V 

+ 9 tt 2 v 2 - 6 tt 2 v 3 - 2 tt 3 v + 3 u 3 v 2 - tt 4 v ) G(l - v, u) H(0, v) i tt 



+ 



39t; + 39f 2 - 13 1; 3 - 13u + 46uv- 53uv 2 + 20 tit; 3 



1 ( 13 _ 

24tt 2 (l-tt)(l-t;) 2 V 

- 13 tt 2 v + 17 tt 2 v 2 - 4 tt 2 v 3 + 6 tt 3 v - 9 tt 3 1> 2 + 3 tt 4 v ) G(l - v, u) H(0, 1, v) 

tt ^-tt \Tt ( 13 - 39 v + 39 w 2 - 13 1; 3 + 5 tt - 8 tt v + u v 2 + 2 tt t; 3 - 9 tt 2 

24 tt 2 (1 - tt) (1 - v) 2 V 

+ 32 tt 2 v - 46 tt 2 v 2 + 23 tt 2 1; 3 - 9 tt 3 + 24 tt 3 v - 18 tt 3 1; 2 + 3 tt 4 v ) G(l - v, u) Q 2 

(277 - 831 v + 831 1; 2 - 277 1; 3 - 223 tt + 850 tt v - 1031 tt t; 2 



+ 



144tt 2 (l-tt)(l-t;) 2 
+ 404tit> 3 - 45 tt 2 - 70 tt 2 v + 140tt 2 t> 2 - 25 tt 2 v 3 + 9 tt 3 + 96 tt 3 v - U4u 3 v 2 - 18 tt 4 
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+ 57u 4 v) G(l-v,u)H(0,v), 



| (3) 



(l-v) 



5 , 11 . 11 rT ,„ w 49 7 , 7 TTln . A 389 A 
+ l^ C2+ l^^ C2+ 16 i/( °' u)C2 -32 C4 -4 i7rC3 -4 i/(0 ' wK3 -144 C3 



2759 2759 „.„ , 81659 

H i 7T H H(0, v) H 

864 864 v ' ' 10368 



4 4) = o, 



(B.2) 



(1) = «(l-u-jj 
7 8u(l-u) 3 



2 - 5 u + 5 



2 i tt C 2 - 2 G(0, u) C2 - G(0, u) i? (1, 0, i>) 



- H(l,v)C 2 - H{l,0,v)in- H(l,l,0,v 
v(l — u — v) 



H - '- G(l —v,l — v, —v, u) i tt — G(l —v,l — v, —v, u) H(0, v) 

+ G(l-v,l- v, -v, 0, u) + G{l-v,l- v, u) ( 2 + G(l - v, 1 - u, u) H(Q, v) i tt 
-G(l-v,l- v, u) H(0, 0, v) + G(l-v,l- v, u) H(l, 0, v) + G(l - v, 1 - v, 0, u) H(0, v) 
+ G(l - v, u) Cs + G(l - v, u) H(0, v) ( 2 + G(l - v, u) H(0, 0,v)in 

- G(l - u, u) ff(0, 0, 0, v) + G(l - u, it) ff(0, 1, 0, v) + G(l - t>, u) i? (1, 0, 0, v) 

- G(l - v, 0, -u, it) i tt + G(l - u, 0, -v, u) H(0, v) - G(l - u, 0, -v, 0, u) 

+ G(l - v, 0, u) £f(0, 0, v) + G(l, 1 - v, -v, u) iir - G(l, 1 - w, -v, u) H(0, v) 

+ G(l, 1 - v, -v, 0, u) + G(l, 1 - v, u) C 2 + G(l, 1 - v, u) H(0, v) i tt 

-G(l,l-v, u) H(0, 0, v) + G(l, l-v,u) H(l, 0, w) + G(l, l-v,0, u) H(0, v) 

+ G(l,u) ( 3 -2 G(l, u) i tt (2 - G(l,u) H(0, v) ( 2 - G(l,u) H(0, 0, v) i tt 

+ G(l,u) H(0,0,0,v) - G(l,u) H(0,l,0,v) - G(l,u) H(l,v) ( 2 - G(l,u) H(l,0,v)iTT 

+ G(l,u) H(l, 0, 0, v) - G(l,u) H(l, 1, 0, v) + G(l, 0, -v, u)i-K- G(l, 0, -v, u) H(0, v) 

+ G(l, 0, -t), 0, w) - 2 G(l, 0, u) C2 - G(l, 0, tt) H(0, 0, v) - G(l, 0, u) i7 (1, 0, u) 

— ( 4 - 10 v + 8 v 2 - 2 v 3 - 12 u + 30 u v - 23 u v 2 + 5 u v 3 + 12 u 2 

Cs 



+ 



8«(1- u) 3 (l- v) 2 

- 34 u 2 v + 26 u 2 v 2 - 5 u 2 v 3 - u 3 + 11 u 3 v - 7 u 3 v 2 - 4 u 4 + u 4 v + u b ^ 

+ G(l - v,-v,u)in- G(l - -u, u) H(0, v) + G(l - u,-i;,0,u) + G(l - v,u)( 2 
+ G(l - u,u) H(0,v)i7T- G(l - i?(0,0,u) + G(l - u,u) ff(l,0,u) 

+ G(l - v, 0, u) H(0, v) + H(1,0, 0, v) 

H 7 ttt To f2w-4w 2 + 2w 3 + 2u-8uw+llMi; 2 -5uw 3 -8u 2 + 16u 2 w 

8tt (1 — m) j (1 — w)^ V 

- 14 u 2 v 2 + 5 u 2 v 3 + 9 u 3 - 9 u 3 v + 3 u 3 w 2 - 4 u 4 + it 4 v + it 5 ^ - G(0, -w, u) i 7r 
+ G(0, -v, u) H(0, v) - G(0, -u, 0, u) + G(0, u) i? (0, 0, v) + i? (0, u) C2 + H(0, 0, u) i tt 
-H(0,0,0,v) + H(0,l,0,v) 

''"'I + 3! ,) g (0, ») + M g _- « - «) (3 - „ - .) f » 



+ 



288(1 - w) 2 



-27- 



+ 



+ 



(1 — u — v) 
48(1 -u) (l-v) V 
1 



16(l-u) 2 (l--u) 2 
5u 2 v 2 -2u 3 -u 3 v^ H(0,v)iir 
1 



(l3 + 6w-19u)G(0,u) 

2 + 3v -7v 2 + 2v 3 -4u-l3uv + 20uv 2 -6uv 3 + 4 u 2 + 7u 2 v 



^6 — 13 v -\- v 2 -\- 6 v 3 — 12 u -\- 27 uv-\- 16 u v 2 — 18 u v 3 -\- 12 u 2 



48 (1 - u) 2 (1 - v) 2 

- 45 u 2 v + 7 u 2 v 2 - 6 u 3 + 19 u 3 v ) G(0, u) i? (0, u) 

) fl44w-144t; 2 -133u+133Mi; + z 

288u(l-u) V 

H 7 To~7 7 f-2i; 2 + 2w 3 -it-uw + 7uu 2 -4uw 3 + 2?i 2 + 2u 2 i;-5u 2 w 2 

8u (1 — u) 2 (1 — v) \ 

- 2 u 3 + u 3 v + u 4 ) - G(-v, u) i 7T + G(-v, u) £f(0, v) - G(-v, 0, u) 

- 771 — 7 \tt, 777 f 6 v 2 - 12 w 3 + 6 v 4 + 3 u - 22 u v - 2 u v 2 + 33 u v 3 - 12 u v 4 

24 m (1 - u) 2 (1 - v) 2 V 

-\- 66 u 2 v — 23 u 2 v 2 — 15 u 2 v 3 + 6 u 3 - 75 it 3 v + 25 it 3 « 2 - 3 u 4 + 25 u 4 v ) #(0, 0, u) 

+ 777 — 7 TT77 77 f 12 v - 36 v 2 + 36 v 3 - 12 v 4 + 6 u - 49 u v + 97 u v 2 - 78 u v 3 + 24 u v 4 

48 u (1 - u) 2 (1 - v) 2 V 

- 12 u 2 + 51 it 2 v - 32 u 2 v 2 + 6 u 2 w 3 + 12 u 3 - 45 u 3 u + 7 u 3 w 2 - 6 u 4 + 19 u 4 v ) 0, i>) 



4 ^2 

6 u 



+ 



24t> - 72 v 2 + 72v 3 - 24v 4 + 6u - 85uv+l93uv 2 - I62uv 3 



48u(l - u) 2 (1 - v) 2 

+ 48uv 4 - I2u 2 + 7Hu 2 v - 80u 2 v 2 + 30u 2 v 3 + 12 u 3 - 45 u 3 v + 7 u 3 v 2 - 6 u 4 + 19 u 4 v ) C2 



Al 2 ) 



w (1 — u — f) 



8li(l-!l) 2 



(2 -3m) [-3C3-*7rC2-G(l-«,«)H(l,0,t;)+4H(l,u)C2 



+ H(l, 0, t>) i 77 -H(l,0,l,v) + H{l,l,0,v) 

v (1 — u — v) 



16u(l -it) 2 



24 it (1- it) 
?; (1 — u — u) 
48 it (1 - it) 2 

l> (1 — U — l>) 



^6— 5u) G(l — u, u) iir — H(l,v) iir 
2 ^22- 39m) ff(l,0,u) 



^26 - 51 it) - G(l - v, 1 - v, it) + G(l - v, u) H(l,v) - H{1, l,v) 

v - u 3 ) 



— (4 - 10 v + 6v 2 - 6u+ lAuv - §uv 2 + 2u 2 - 



\u{\- u) 2 (1 - v) 

G(l - u, 1, 1 - v, u) + G(l - 1, u) i 7T - G(l - w, 1, u) H(0, v) - G(l - u, 1, u) i? (1, v) 

v (1 — u — v) 



— - ( 2 v - 2 v 2 - 2 u v + 3 u v 2 - 2 u 2 v + u 3 
v) z \ 
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-H(0,0,v) 



u(l-u) 2 (l-v) 

+ G(0, l-v,u) H(0, v) + G(0, l,l-v,u) + G(0, l,u)in- G(0, 1, u) £T(0, u) 
-G(0,l,u)H(l,v) 
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+ 



+ 



v (1 — u — v) 



48u(l-u) 2 (l-u) 



^—2 ^44- U4v + 70v 2 -66u+158uv- 105 u v 2 + 26 u 2 v - 13 u 3 ^j 



H(0,l,v) 



v (1 — u — V 



-v) 2 V 



44 - 70 w + 26 v 2 - 66 u+ 114 wu- 39 uw 2 - 18w 2 v 



48u(1-u) 2 (1-i;) 2 
+ 9« 3 ) if(0,u)i7r 

~ ,n 7 ^ ~ "oT-,^ \ 9 (l60 - 294 v + 134 « 2 - 228u + 430uw- 189 uw 2 - 26 u 2 w + 13 u 3 ) ( 2 
48u(l - u) 2 (1 - v) 2 V / " 



+ 



w (1 — u — v) 
4u 2 



y G(l, «) if (0, v) i 7T - y G(l, u) if (0, 1, v) + y G(l - u) if (0, 0, v) 
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- — G(l, u) if (0, 0, v) + 2 G(l - w, 1 - v, u) C2 + 2 G(l - v, 1 - v, u) if (0, v) i ir 

- 2 G(l - v, 1 - v, u) if(0, 1, v) - 2 G(l - v, 1 - v, 1, 1 - u, u) - 2 G(l - v, 1 - v, 1, u) i tt 

+ 2G(l-t),l-«,l,u) if (0, !))+2G(l-t),l- w, 1, u) if (1, v) + 2 G(l - v, 0, 1 - u, u) if (0, u) 
+ 2 G(l - w, 0, 1, 1 - i), u) + 2 G(l - v, 0, 1, u) i tt - 2 G(l - v, 0, 1, u) if (0, v) 

- 2 G(l - w, 0, 1, u) if (1, u) + G(l - v, 1, 1 - u) if (0, u) - 4 G(l - v, 1, u) C2 

- G(l - u, l,u)if(0,u)i7r + G(l - l,u) if(0, l,u) - G(l - w, l,u) H(l,0,v) 
+ 2G(l-v,l,l,l-v,u) + 2G(l-v,l,l,u)in-2G(l-v,l,l,u)H(0,v) 

- 2 G(l - u, 1, 1, u) if (1, u) + 2 G(l, 1 - v, 1 - w, u) if (0, v) - 4 G(l, 1 - u, m) C2 

- G(l, 1 - u, u) if (0, u) i tt + G(l, 1 - v, u) if (0, 1, v) - G(l, 1 - v, u) if (1, 0, u) 
+ 3 G(l, 1 - v, 1, 1 - i>, u) + 3 G(l, 1 - v, 1, u) i tt - 3 G(l, 1 - v, 1, u) ff(0, u) 

-3G(l,l-u,l,«)H(l,i;)-3G(l,u)C3-G(l,«)*7rC2+4G(l,«)H(l,u)C2 
+ G(l, u) H (1, 0, u) i tt - G(l, u) ff(l, 0, 1, v) + G(l, u) ff(l, 1, 0, u) 

- G(1,0,1- u,u)if(0,u) - G(1,0,1,1-u,m) - G(1,0, ztt + G(l, 0, 1, u) if(0, v) 
+ G(l,0,l,u)if(l,u) + 2G(l,l,l-u,l- !),u)+2G(l,l,l- v,u)i7r 
-2G(l,l,l-w,u)if(0,w)-2G(l,l,l-u,u)if(l,w)-2G(l,l,u)C2 

- 2 G(l, 1, u) if(l, u) i tt + 2 G(l, 1, u) ff(l, 0, v) + 2 G(l, 1, u) ff(l, 1, u) 

-2G(1,1,1,1- -2 G(l, 1,1, w)i7r + 2 G(l, 1,1, u)H(0,v) + 2 G(l, 1,1, u) if (l,u) 

+ ^iT" - ^ f 3 - 2 u - 3 u 2 ) [G(l, 1 - u, u) i tt - G(l, u) H(l,v) i tt 
8u z (1 — u) z \ / L 



v (1 — m — v) 



12 u 2 (1 - 

v (1 — 11 — 1 



^11 -28u + 20?i 2 ) G(l, 1 - v,u) H(0,v) - G(l,u) if (1,0, w) 



^13 - 38 u + 31 ?i 2 ) - G(l - 1 - v, u) if (0, d) - G(l, 1 - v, 1 - u) 



24u 2 (l-u) 2 
+ G(l, l-v,u) H(l,v) - G(l,u) ff(l, l,u) 
t; (1 — u — v) 



+ 



24u 2 (l-u) 2 

+ G(l,l,u)if(l,w) 

v (1 — u — v) 
~ 24u 2 (l-w) 2 



^13 - 14 u - 5 u 2 ) - G(l, 1, 1 - v, u) - G(l, 1, u) i w + G(l, 1, u) if (0, u) 



^67- 122?i + 49u 2 ) G(1,m)C 2 
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+ 7, — tt~~ — !t ; ^ „ ( 3 - 6 v + 3 v 2 - 6 u + 14 u v - 8 u v 2 + 3 u 2 - 8 u 2 v + 6 u 2 v 2 
8 u (1 — u) 2 (1 - v) 2 V 

- 2 u 3 v + u 4 ) G(l - w, u) if(0, v) i tt 

v(l-U-v) f 10 ofi ^ | 10 ^ 2 „«„,„, on„.„,2 



+ 



24u 2 (l-u) 2 {l-v) 



-v) 2 ( 



13 - 26 v + 13 v l - 26 u + 46 u w - 20 u v l + 13 h 2 - 20 u 2 v 



+ Au 2 v 2 + 6u 3 v-3u 4 ^jG(l-v,u)H(0,l,v) 
v (1 — u — v) 



+ 
+ 
+ 

+ 
+ 



. . 13 - 26 v + 13 v 2 - 8 u + 10 u v - 2 u v 2 - 14 u 2 + 34 u z v 
24 u 2 (1 - u) 2 (1 - w) 2 V 

23 u 2 v 2 + 6 u 3 v - 3 u 4 ) G(l - w, u) C2 

— ^7-— ——^-77 ( - 554w + 554 + 36u + 987uv - 909uv 2 - 150u 2 - 78u 2 v 

288 u (1 — u) 2 (1 — f J 



114 u- 



G(l - + H(l,w) 



(1 — ti — f ) 



288u(l-u) 2 (l--y) 
(1 — u — v) 



f554u 

- v) V 



554 w 2 - 36 u - 855 u v + 909 u v 2 + 18 u 2 - 54 u 2 w + 18 w : 



— (554v 2 - 554v 3 + 96uv- 842uw 2 + 909uv 3 -132u 2 + 228u 2 t> 
288u(l-it) 2 (l--y) 2 V 

422 u 2 v 2 + 132 u 3 + 31 ?i 3 v ) i? (0, u) 



(1 — u — v) 



277 v - 277 v 2 -18u- 6Uuv + 614ui; 2 + 36u 2 + 358 u 2 v 



144w 2 (l -u) 2 (l-v) 
340 u 2 v 2 - 36 u 3 - 18 u 3 v + 18 u 4 ) G(l, 1 - w, u) + G(l, it) i tt - G(l, it) ff(0, w) 



+ 



—-r(- 277v + 554 w 2 - 277 v 3 + 18 u + 464 u v - 886 uv 2 + 404 u w 3 
f 2 V 



+ 



G(l,«)H(l,t;) 

(1 — u — v) 
144u 2 (l -w) 2 (l -v) 

36 u 2 - 124 u 2 v + 224 u 2 v 2 - 25 u 2 v 3 + 36 it 3 - 18 u 3 v - 96 u 3 v 2 - 18 u 4 + 57 u 4 v ) 
-G(l - w,u)i?(0,-y) 

(lUv - 144 w 2 + 144 u- Uuv- 133uv 2 -421u 2 + 144 m 2 w + 277 u 3 ) , 

u 1 — v V / 



288«(l-u)(l-t;) 

4 3) = o, 



4 4) = o. 



(B.3) 



C. Helicity amplitudes for qq — > Vg and — > Vq 

In [21] the helicity amplitudes for the processes 

7*(j> 4 ) — > q(pi) + q{p2) +g{p3) (C.l) 

have been computed. As already discussed, the kinematical region relevant for 3-jet pro- 
duction is characterised by q 2 and Sij all positive. 
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Following [24], we can analytically continue these matrix elements to the kinemati- 
cal configuration relevant for (V + lj) production at hadron colliders, which is the same 
kinematical situation relevant for the production studied above: 

?(P2) + g(Pi) — »• g{-Pa) + l*{Pi) , (C2) 

and 

q{p2) + g(p 3 ) — »• g(-Pi) + 7>4) , (C.3) 

The kinematical situation in (C.2) is the same as for, 

g(P2) + ?(Pi) — »• 7(-Ps) + V{ PA ) , (C.4) 
so we can introduce again the dimensionless variables 

»13 y q 2 I 

u = = — , v = — = — , 

S\2 X S 12 X 

which fulfil 

0<u<v, < v < 1 . 

On the other hand, the kinematical situation in (C.3) can be described with the following 
choice of dimensionless variables [24]: 

/ «i3 y / q 2 i fr ,^ 

u = = — , v = — = - , C.6 

which fulfil again: 

< u < v' , < v' < 1 . 

The helicity amplitude coefficients a, j3 and 7 are vectors in colour space and have 
perturbative expansions: 



n = v / 47ra v / 4vra s Tfj 



for Q = a, (3, 7. The dependence on (it, v) or (1/, 1/) is again implicit. 

The ultraviolet and infrared properties of the helicity coefficients are fully described 
in [21], 

n (0) = j 2 (0),un ) 



(C7) 



2e 

n (2) = Se -2 n (2),un _ 3^ Sr l n (l),un _ ^ _ n (0),un 

and 

= /«(e)n( ' + nW' finite , 

f}( 2 ) = Ll/W (e )jW (e ) - ^jW (e ) + e-^ ^" 2 ^ + K) /«(&) + H^(e)\ n<°> 
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+I^\e)Q^ +Q^ Mite , (C.8) 

where the infrared singularity operator I^(e) is a 1 x 1 matrix in colour space and is given 
by 



2T(1 - e) ^ 
where (since we have set /U 2 = S123 



(C.9) 



s« = (-—)• (CIO) 



Expanding Sij, imaginary parts are generated, the sign of which is fixed by the small 
imaginary part +i0 of Sij. 
Finally we have 



" ,2, < £ > = 47Tf^^ ,2, • < c -"> 



where, in the MS scheme: 



V 3 432 72 J V 2^ 3 54 48 7 V 16 4 J N 2 

V 18 36/ V 54 24/ TV 27 F v ' 

They finite one-loop amplitudes can be decomposed according to their colour structure 
as follows: 

nW' & ™ te (u,v) = Nan(u,v) + -bn(u,v) + p cn(u,v) . (C.13) 

We attach to the arXiv submission of the paper the one loop coefficients expanded up to 
0(e 2 ). 

The finite two-loop remainder is obtained by subtracting the predicted infrared struc- 
ture (expanded through to C(e )) from the renormalised helicity coefficient. We further 
decompose the finite remainder according to the colour structure, as follows: 

n (2),finite( U)t; ) = N ^A n (u,v) + B n {u,v) + -Lc n (u,v) + NN F D n (u,v) 
Nf„ , , Ar2 ^ , „ ( 4 



+ -Jj-En(u,v) + N F F Q (u,v) + iV F , y I — - iVj G n («,v) , 

The complete expressions in FORM format can be found attached to the arXiv submission 
of this paper. 
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